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Meson widths from string worldsheet instantons 
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We show that open strings living on a D-brane which lies outside an AdS black hole can tunnel 
into the black hole through worldsheet instantons. These instantons have a simple interpretation 
in terms of thermal quarks in the dual Yang-Mills (YM) theory. As an application we calculate the 
width of a meson in a strongly coupled quark-gluon plasma which is described holographically as a 
massless mode on a D7 brane in AdS^ x S5. While the width of the meson is zero to all orders in the 
1/VA expansion with A the 't Hooft coupling, it receives non-perturbative contributions in l/vA 
from worldsheet instantons. We find that the width increases quadratically with momentum at large 
momentum and comment on potential phenomenological implications of this enhancement for heavy 
ion collisions. We also comment on how this non-perturbative effect has important consequences 
for the phase structure of the YM theory obtained in the classical gravity limit. 



A heavy quarkonium bound state, Hke J /ip or T, which 
finds itself in the quark-gluon plasma (QGP), becomes 
increasingly unstable and eventually dissociates at suffi- 
ciently high temperatures. On the one hand this can be 
attributed to the weakening attraction between a heavy 
quark and anti-quark in the bound state due to color 
screening of the medium [1]. On the other hand the 
bound state can be broken up from collisions with the 
deconfined quarks and gluons in the medium [2]. Given 
that the QGP at RHIC and LHC is likely strongly cou- 
pled, understanding such medium effects on the propa- 
gation and dissociation of heavy quarkonia presents non- 
trivial theoretical challenges which must be confronted 
in order to interpret experimental data on quarkonium 
suppression. 

Interesting insights have recently been made into this 
problem from strongly coupled model gauge theories 
using the AdS/CFT correspondence. AdS/CFT-bascd 
methods are powerful at attacking questions of dynami- 
cal origin, such as how the motion of quarkonia relative 
to the medium affects their various properties. The sim- 
plest example of the correspondence is provided by the 
duality between TV = 4 SU{Nc) super Yang-Mills (SYM) 
theory and string theory in AdS^ x 5*5 [3] . Based on a cal- 
culation of the potential between a heavy external quark 
and antiquark pair moving in the strongly coupled hot 
N — 4: plasma, it has been argued in [4] (see also [5] ) that 
the dissociation temperature Tdiss of a heavy quarkonium 
decreases with their velocity v relative to the medium as 
Tdiss{v) Ki {1 — v^)iTdiss{v — 0). Such a velocity scaling, 
which can be heuristically understood as due to increased 
screening from the boosted medium, could provide a sig- 
nificant additional source of quarkonium suppression at 
nonzero transverse momentum in heavy ion collisions [4] . 

Rather than drawing inferences from the heavy quark 
potential, it is also possible to directly study the prop- 
agation of mesons in a strongly coupled plasma. While 
A/" = 4 SYM theory itself does not contain dynamical 
mesons, one can obtain a closely related theory which 
does contain mesons by adding to it Nf <^ Nc funda- 
mental "quarks" , which corresponds to adding some D7- 



branes to AdS^ x 6*5 in the gravity picture [6]. It was 
found in [7, 8] that meson dispersion relations are dra- 
matically modified by the plasma and in particular, there 
exists a limiting velocity Vc{T) < 1, which decreases with 
increasing temperature. The existence of a subluminal 
limiting velocity is consistent with the velocity-enhanced 
screening obtained from the heavy quark potential, as 
when V > Vc{T) the quark and anti-quark are completely 
screened and no bound states are possible. 

For a more complete understanding of the dissocia- 
tion of mesons one needs to study their widths. We 
will be particularly interested in the momentum depen- 
dence of the widths. This has not been possible within 
the classical gravity approximation developed so far. In 
this approximation, the mesons are stable (i.e. they have 
zero width) for T smaller than a dissociation temperature 
Tdiss, but completely disappear for T > T^iss [7, 9, 10]. 
The approximation also has another important draw- 
back: the densities of quarks and antiquarks are iden- 
tically zero for a range of temperatures and chemical po- 
tentials [11] even though they should obey the standard 
thermal distribution. 

In this paper, we discuss a novel tunneling effect on 
the string worldsheet which gives rise to nonzero quark 
densities and meson widths for T < Tdiss- This enables 
us to calculate explicitly the momentum dependence of 
the width of a meson in a strongly coupled QGP. We find 
that the width increases quadratically with momentum 
at large momentum. 

At finite temperature. A/" = 4 SYM theory can be de- 
scribed by a string theory in the spacetime of a black 
hole in AdSs x 6*5, whose metric can be written as 
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where f = ^ [1 — ^] , x = (xi, 0:2, X3). dfl^ is the met- 
ric on a unit five-sphere 6*5 which can be written as 
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sin^ edill 
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(2) 



with dri| the metric for a three-sphere 5*3. The string 



Note that A^ is a dimensionless number of 0(A°) deter- 



^T^gvM and the curvature radius R is related to the 



9ymNc by -^ = \/A. The per- 



couphng Qs is related to the Yang-Mills^coupling gyu by mined by the ratio L/tq, and m\ ' can be interpreted as 

the effective mass of a quark at temperature T . 

The mesons corresponding to massless fluctuations on 
the D7-brane can be found by solving the linearized equa- 
tions resulting from expanding the DBI action around 
the embedding. For example, the quadratic action for 
the fluctuation x'^(C") of the location of D7-brane in the 
6 direction can be written as 



JYM 

't Hooft coupling A 

turbative gg and a' expansions in the bulk string theory 
are related to the 1/iVc and -j= expansions in the Yang- 
Mills theory respectively. The temperature T of the YM 
theory is given by the Hawking temperature of the black 
hole, T = ■:^^- Adding Nj fundamental "quarks" can 
be described in the dual string theory by adding Nf D7- 
branes in (1) [6]. A fundamental "quark" in the YM 
theory can be described by an open string with one end 
on the D7-branes and the other end on the black hole. 
Open strings with both ends on the D7-brancs can be 
considered as "bound states" of a quark and antiquark, 
thus describing meson-type excitations in the YM theory. 




FIG. 1: An embedding of the D7 brane (green) in the 
AdS^ X Ss black hole geometry for T < Tdiss which lies en- 
tirely outside the black hole. Inset: the Euclidean r — r plane 
aiO = Q showing a world-sheet instanton (red) connecting the 
tip of D7 brane r = r^ to the horizon at the center of the 
disk r = ro. 

We now briefly outline the standard procedure for ob- 
taining the meson spectrum [9]. We will take Nf = 1, 
Nc — > oo, and A large but finite throughout the paper. 
The D7-brane can be chosen to lie along the directions 
^" = {t,x,fl3,9) and using the symmetries of the prob- 
lem the embedding in the two remaining transverse di- 
rections can be taken as 0(^") = and r(^") = r{9). 
At the lowest order in the a' expansion, r{d) can be 
determined by extremizing the Dirac-Born-Infeld (DBI) 
action of the D7-brane with the boundary condition 
r{9) cos 6\g^iL — > L, where L is related to the mass rriq 
of a quark in the Yang- Mills theory as niq ~ 2^^- For 
T smaller than some Tdiss, fi^) has the form shown in 
Fig. 1. The brane is closest to the black hole at 6* = 0, 
where there lies a 4-dimensional subspacc spanned by 
(t, x) since here the 5*3 in (2) shrinks to a point. De- 
noting rm = '''(O = 0) > ro, the shortest open string 
connecting the D7-brane to the horizon has a mass in 
the YM theory 
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5,,B7[x1 = -f / A 



where /iy 
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gG^^g''^daX%X^ (4) 
is the tension of the D7-brane, 



Gflxf, = R cos 9, and g^p denotes the induced metric on 
the D7-brane. Writing x"^ = e-*'^*+'^-^yi(r23)V'(6'), the 
equation of motion for "0 can be written as 



H{k,l)ij{9) =uj'^il^{e) 



(5) 



where H{k,l) is a second order differential operator in 
and Yi are spherical harmonics on the 6*3. For a given 
fc, I, H{k, I) has only discrete eigenvalues w^ labeled by an 
integer n, giving rise to dispersion relations uj = cij„(fc, /), 
all of which have zero width. In particular, the meson 



masses are of order AI = 



_ 2V2L _ 4\/27r 



Since M is 



parametrically smaller than niq in vA, the mesons have 
a large binding energy, given by 2mq . There exists a 
temperature Tdiss = 0.122M, beyond which the D7 brane 
falls into the black hole and mesons cease to exist as well- 
defined quasi-particlcs [7, 10]. 

We stress that the zero-width conclusion only depends 
on the topology of the embedding in Fig. 1. Since mesons 
can only dissociate by falling into the black hole, when 
the D7-branc lies above the black hole horizon the mesons 
are necessarily stable. Given that the brane embed- 
ding and the background geometry are smooth, includ- 
ing higher order perturbative corrections in a' should not 
change the topology of the brane embedding if the dis- 
tance between the brane and the horizon is parametri- 
cally larger than the string scale. This implies that the 
widths of mesons should remain zero to all orders in the 
perturbative -j= expansion. 

One can also turn on a quark chemical potential /i < 

(T) 

niq in the boundary theory by setting At = fi, where 
At is the time component of the gauge field on the D7- 
brane [11, 12]. Since the DBI action and its higher or- 
der a' corrections contain only derivatives of At, the D7- 
embedding and the meson spectrum are not modified by 
turning on the constant mode of At. Thus, the meson 
widths and the net quark density remain zero to all or- 
ders in the a' expansion even at a finite chemical poten- 
tial [18]. 

The above conclusions can be further illuminated by 



simple thermodynamic reasoning. From (3), (3m, 



(T) 



oc 



vA, the quark (or anti-quark) density, being proportional 
to e"'^™'' '^'^j is then exponentially suppressed in vA 

(T) 

for /i < rriq . Similarly, since the binding energy Ebe of 

(T) 

a meson is 2mq , thermal effects which destabilize the 
mesons are also exponentially suppressed in vA. Thus 
the meson widths and quark densities are not visible in 
the perturbative expansion in 1/vA and can only have 
non-perturbative origins on the worldsheet. 

There are indeed non-perturbative corrections in a' 
which effectively change the topology of the D7-brane 
embedding and generate non-vanishing meson widths 
and quark densities. To see this it is more convenient 
to analytically continue (f ) to Euclidean signature with 
t -^ —IT. Then the r — t plane has the topology of a 
disk. The angular direction r has a period given by the 
inverse temperature (3. The center of the disk is located 
at r = rg. Open strings on the D7-brane are described 
by worldsheets with the topology of a disk whose bound- 
ary lies on the D7-brane. Denoting the worldsheet co- 
ordinates as p € [0,1] and a G [0,27r), the worldsheets 
separate into different topological sectors corresponding 
to the winding number m of the target space disk (r, t) 
wrapping around the worldsheet disk {p,(y). The DBf 
action arises from the sector of trivial winding number 
771 = 0, in which [p, a) is mapped to a single point on 
the D7-brane. In all the other (nontrivial) sectors, the 
string worldsheet is mapped to the region in the r — t 
plane from the location of the D7-brane all the way to 
the horizon r = vq (see inset of Fig. 1). When analyti- 
cally continued back to the Lorentzian signature, such a 
worldsheet describes a tunneling process in which a tiny 
neck is generated between the D7-brane and the black 
hole horizon. As a result mesons can leak through the 
tiny neck into the black hole and dissociate. 

As an illustration, we now calculate the contributions 
from 771 = ±1 sectors to the quark density and the widths 
of mesons in (4). We will only be interested in the lowest 
order term in the a' expansion. The relevant spacetime 
effective action for the D7-brane can be obtained from 
the worldsheet path integral [15] 



SeIx'^ 



DXe 



-nX] + f^^,dat.j^-lB[x^] 



(6) 



disk 



where X = (^q, , r, (/)) denotes collectively the worldsheet 
fields. I[X] is the worldsheet action, which for our pur- 
pose can be taken to be the Nambu-Goto action for a 
string propagating in (1) 



I[X] = J / dadp ^/dethab 



(7) 



with hab = GMNdaX^'dbX" the induced metric on the 
worldsheet and Gmn the Euclidean version of the metric 
(1). The second term in the exponential of (6) corre- 
sponds to turning on A^ = —ip which gives in the Eu- 
clidean theory a nonzero (quark) chemical potential p in 



the boundary YM theory. IbIx"^] = §p^ida ^x^dp^ 
is the boundary action which couples the worldsheet to 
X'^(^"). We have suppressed any dependence on space- 
time and world sheet fermions. The integral (6) can 
be evaluated using the saddle point approximation in 
each topological sector [16], i.e. Se — Sm=o + Sm=+i + 
Sm=-i + ■ ■ ■ with Sm=0 = Sdbi- 

For 777 = ±1, (7) has a classical solution given by 






r^f{p), = 0, 



0, X ~ x^ (8) 



where xq is an arbitrary constant position vector and 
f{p) is chosen so that f(l) = r,„ and /dr^ -I- ^df^ oc 
dp^ + p^da^. Eq. (8) represents the worldsheet of a string 
connecting the tip of the brane to the horizon with the 
± sign corresponding to opposite orientations. It has a 
classical action I± = (3mq ' where 777, ' is the effective 
quark mass introduced in (3). One can readily verify that 
(8) minimizes the action and satisfies the right boundary 
conditions at the D7-brane. Note that there are only 
three bosonic zero modes in (8) [19], since it costs energy 
to move away from 9 = and the worldsheet time a now 
coincides with r. With x"^ set to zero, Eq. (6) then yields 



^TTl — 



m=±l 



e-/3m<-)g 



9s 



(9) 



where the e^^ arises from the second term in the ex- 
ponential of (6), V3 is the spatial volume from integrat- 
ing over the three zero modes in (8), and the — fac- 
tor arises because we are evaluating the disk path inte- 
gral. Z? is a real number coming from Gaussian integra- 
tion around the saddle point (8) (including worldsheet 
fermions) whose sign we will fix from physical require- 
ments. Identifying the Euclidean action with (3F{(3, p) 
where F{(3,p) is the free energy, equation (9) leads to a 

net quark charge density — ^e~'^"'<! sinh/3/i, which in 
turn requires that D should be negative [20] . It is natural 
to interpret 5'm=±i as the contributions from quarks and 
anti-quarks separately: Sm=±i = —n±V3, which from (9) 
leads to a quark and antiquark number density given by 
77± = e-'^™f'e±^'3_L (_£)), Note that 77± ex l/g., oc N^ 
since quarks come in an A^c-niultiplet. 

In our derivation of (9) we have assumed the embed- 
ding of the D7-brane is given by that determined by the 
DBI action. This is justified for p < nVq ' since the cor- 
rection to the DBI action is exponentially small. When 

(T) 

P > "m-q , the backreactions from instantons become 
large and the embedding of Fig. 1 cannot be trusted. 

The nonzero quark densities for nonzero p < rriq have 
important implications for the phase structure of the the- 
ory. As discussed in [11] (see also [12]) based on the anal- 
ysis of the DBI action (which corresponds to A = oo), at 
low temperature there is a phase transition in chemical 
potential at which the net quark charge density becomes 



nonzero. Our results strongly indieate this phase transi- 
tion (at nonzero temperature) is smoothed to a crossover 
at any finite value of A. 

To find the widths of the mesons described by (4), 
we need to compute (6) to quadratic order in x'^ ■ For 
simplicity, we will restrict to the / = mode on the 
S'3. Near 6* = 0, the worldsheet action for (j) is given 
by -^^ J dpda (90)^, which is free. The path integral is 
then straightforward to compute and yields for Sm=±i 



The ratio 



d^xodrdr' (xHr,So)GDi<y,<j')x*{r',So) 



R^n± 

27ra'P'^ J " \' " ' - ■ "-/6(=o 

(10) 
where Gd{o', <y') — limp_»i p'_,i dpOpiGoip, cr; /?', a') with 
Gd{p, c; p' , a') the Dirichlet propagator for a canonically 
normalized massless field on the unit disk and a = ^t. 

Note that (10) only depends on the value of x'^ at the 
tip of the brane and is nonlocal in the Euclidean time 
direction. 

Treating (10) as a small perturbation to (4), one can 
compute the corrections to the Euclidean two-point func- 
tion of the (meson) operator dual to x'^ in the boundary 
YM theory, from which the corrections to the Lorentzian 
retarded function can be found by analytic continuation. 
One finds that the poles of the retarded function now 
obtain a nonzero imaginary part. Alternatively one can 
directly obtain the Lorentzian counterpart of (10) by an- 
alytically continuing the worldsheet disk to Lorentzian 
signature with a = irj = i27rt//3, which gives the part 
of Rindler spacetime ds^ = dp^ — p^drf with p < 1. 
The Lorentzian spacetime effective action can be ob- 
tained using the Schwinger-Keldysh contour, giving the 
Lorentzian equation of motion [21] 



do. {V^gG^^d^x^) 



Rh-i±s{e) 



dt' GR{Ti^r/)xHt') 



where Gr^ii — ti') = limp_»i,p' 



(11) 
1 dpdp'GR{p,r];p',r]') with 
Gr{p, 77; p' , T]') the retarded propagator for a massless field 
in the Rindler spacetime with Dirichlet boundary condi- 
tion at p = 1. Fourier transforming (11) to momentum 
space and using j drie^^^GR{ri) = w, one finds that (5) 
is modified to 



rn(fc) |V„(0 = O;fc)|^ 



rn(0) |V„(^ = 0;fc = 0)|2 



(14) 



can be evaluated numerically and the results are shown 
in Fig. 2. For large k, the asymptotic form of the wave 
function, found in [8], can be used to show that the 
width (14) scales hke fc^ for large k: r„(fc)/r„(0) = 
Rn[T/M][k/Mf + 0[k) for some function Rn[T/M]. 
Furthermore for temperatures T -^ M and k ^ 



M-' 



one finds the closed form expression 

2(47r) 



r.(o) 



T-fc^ 



7 — rrrr? — „ ;n\ . rfi ■ Fig. 2 has the interesting feature that 
the width is roughly constant for small fc, but turns up 
quadratically around k/M = 0.52(Tdiss/T)^, which is 
roughly the momentum at which the meson approaches 
its limiting velocity Vc{T). This is consistent with the 
conclusions based on the velocity dependence of the 
screened quark potential found in [4]. Note that the 
width as defined here is in the rest frame of the plasma, so 
the fc^ behavior at large k should be contrasted with the 
1/fc behavior of a relativistic decay width which comes 
from the usual time-dilation argument. 

The plots here also share some similarities with those 
in [13] for momentum-dependence of meson widths ob- 

(T) 

tained for jj, > rriq with A = 00 where the relevant D7 
brane embedding resembles a long spike reaching down 
to the horizon. 




2 3 

Log(l+ k/M) 

FIG. 2: The behavior of the width as a function of k for 
T/Tdiss = .99, .71, .3, .13 from left to right. The dashed lines 
are analytic results for large momenta. 



H{kJ = 0)^- 






5{e)ip{e)=uj^il;{d) (12) 



with A = ^— ,g(— g"). Writing the dispersion relation 
as cij = Lu„ — 4r„ where n denotes the excitation num- 
ber, and using first order perturbation theory in n± we 
find [22] 



ri-) = '-^^IMO 



NclTl^ 



■n± 



(13) 



with il}n{& = 0) eigenfunctions of (5) evaluated at the tip 
of the brane. Recall that n± are thermal densities for 
quarks and antiquarks and are proportional to Nc- 





FIG. 3: Schematic diagrams of the relevant thermal processes 
contributing to the meson width. For large A the first process 
is dominant, coming from the single instanton sector. 

Our result (13) has a very simple physical interpreta- 
tion as shown in the left plot of Fig. 3, in which a meson 



breaks apart by colliding with a quark (or anti-quark) in 
the thermal bath. There are also processes correspond- 
ing to a meson breaking up by colliding with gluons in 
the thermal bath, shown in the right plot of Fig. 3. For 
such a process to happen the gluon should have an en- 
ergy above the binding energy of the meson. The density 
of such gluons is thus suppressed by e~'^^™i and should 
be described by an instanton and anti-instanton so that 
the resulting worldshcct has trivial topology. Wc expect 
that contributions from such processes are also controlled 
by the value of the meson wave function at the tip of the 
brane, and possibly have similar growth with momentum. 

Our method should be generic to any theory with a 
holographic dual. While the precise value of the width is 
clearly model dependent, it is conceivable that the mo- 
mentum dependence may apply in a wider context in- 
cluding QCD. In particular, our result highlights the con- 
tributions to quarkonium suppression from the collisions 
with medium quarks and gluons in the large transverse 
momentum region in heavy ion collisions. 

Consider the effect of such a momentum scaling on 
J/i/j with AI « 3GeV. The dissociation temperature 
from the gravity set-up is Tdiss = 0.122M « 370MeV 
in fairly good agreement with lattice data [17] for QCD 
Tdiss ~ 2.1Tc for Tc = 170MeV [11]. If we take the 
RHIC temperature of T = 250MeV (this corresponds to 
the second curve from the left in Fig. 2) then a moving 
J/^'s width will increase by a factor of 2(10) at a mo- 
mentum k = 6(13)GeV. When the width of a meson 
approaches the spacing between different meson states, 
the meson can no longer be considered as a well-defined 
quasi-particle. The momentum scaling thus implies a 
maximal momentum beyond which the meson no longer 
exists. As an illustration, suppose the width for the J/ip 
in the QGP at zero momentum is about 200 MeV (which 
is not known) then one expects a maximal momentum 
around 7 GeV. 

Finally, we expect the worldsheet instantons found 
here to have many other applications to various aspects 
of flavor physics in AdS/CFT. 
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Above n > rriq a new phase, where the D7 brane falls 
into the horizon, is thermodynamically preferred. This 
phase exhibits both meson widths and finite quark den- 
sity [11-14]. 

In contrast, the n = sector has eight zero modes cor- 
responding to all directions on the D7-brane. There are 
also no fermionic zero modes here. 

There are a few other indications that D should be nega- 
tive. The instanton action also induces a tadpole for the 
r component of the transverse fluctuations. One flnds for 
D negative the tadpole pulls the brane toward the hori- 
zon as required by physical consistency. Also only for D 



negative, do the meson widths we calculated below have 
the correct sign. 
[21] Eq. (11) can also be obtained directly from the Euclidean 
action (10) using the following general prescription: write 
down the equation of motion following from the Eu- 
clidean action; replace the Euclidean worldsheet time by 
the Lorentzian worldsheet time and the Euclidean prop- 
agator by the corresponding retarded propagator in the 



Lorentzian signature. 
[22] We normalize the eigenfunctions ^„ of H{k,l) as 
-rh^ !^ de,f^G4,4,{-g'')\il}„{k,l-e)\^ = l so that i/;„ 
is dimensionless and has a smooth zero temperature limit. 
It depends on ratios of T, M and k but not explicitly on 
A. 



